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O ' Abstract 



The viruses studied are genetically engineered, charged, semiflexible filamentous bacteriophages 

that are structurally identical to M13 virus, but differ either in contour length or surface charge. 

^ '■ 
r^ , While varying contour length (L) we assume the persistence length (P) remains constant, and 

o ■ 

CJ ' thus we alter the rod flexibility (L/P). Surface charge is altered both by changing solution pH 



and by comparing two viruses, fd and M13, which differ only by the substitution of one charged 
for one neutral amino acid per virus coat protein. We measure both the isotropic and cholesteric 
coexistence concentrations as well as the nematic order parameter after unwinding the cholesteric 



f^ ■ phase in a magnetic field. The isotropic-cholesteric transition experimental results agree semi- 

^^ ■ quantitatively with theoretical predictions for semiflexible, charged rods. 



PACS numbers: 64.70.Md, 61.30.St 



I. INTRODUCTION 

For a suspension of rigid rodlike particles, Onsager determined that hard-core interactions 
alone are sufficient for inducing an entropy driven phase transition from an isotropic phase, 
in which the particles are randomly oriented, to a nematic phase, in which the orientation 
of the particles is distributed about a preferred directionp|. When the rodlike particles are 
semifiexible and/or charged, like many biopolymers such as DNA and F-Actin, the properties 
of the phase transition can differ significantly from thosepredicted for hard, rigid rods. Small 
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amounts of flexibility are predicted |^ and observed [3|, \^ to increase the stability of the 
isotropic phase and lead to a less ordered nematic phase. In this paper, we study the effects of 
flexibility on the isotropic-nematic (I-N) transition using suspensions of the rodlike charged, 
semifiexible M13 virus and M13 virus length-mutants. By varying the contour length (L) 
of our experimental charged rods while maintaining a constant persistence length (P) we 
change the rod flexibility {L/P). The persistence length is defined as the length over which 
tangent vectors along a polymer are correlated JS]. In our experiments the flexibility of the 
rods remains within the semifiexible limit, where P ~ L. The effect of surface charge on 
the I-N transition of charged rods is also investigated. Surface charge is varied by modifying 
both the surface chemistry of the rods and the solution chemistry, by changing pH. 

While Onsager developed the original theory for the isotropic-nematic transition of hard 
and charged rigid rodlike particles, Khokhlov and Semenov were responsible for incorpo- 
rating flexibility into this theory j^. They extended Onsager's theory to include systems 
of semifiexible rods with a large length (L) to diameter (D) aspect ratio (L/D) and arbi- 
trary persistence length. They explicitly calculated the equilibrium properties of the I-N 
phase transition in the limit of very flexible L/P » 1 and very rigid L/P << 1 rods and 
interpolated between the two limits to find the properties of semifiexible rod phase behav- 
ior. Shortly afterwards, Chen numerically calculated the concentrations of the coexisting 
isotropic and nematic phases as well as the order parameter of the coexisting nematic phase 
for arbitrary flexibility using Khokhlov- Semenov theory 6|. For rigid rods, the limit of sta- 
bility of the isotropic phase is predicted to be q = 4/6, where q is the number density 
and b = tcL'^D/A, the average excluded volume in the isotropic phase |2|. For flexible rods, 



Khokhlov-Semenov theory predicts that shght semiflexibihty will increase the stability of 
the isotropic phase by increasing bci, and will narrow the I-N coexistence region. Flexibility 
is also predicted to significantly lower the nematic order parameter at coexistence. The ne- 
matic order parameter S is the second moment of the orientational distribution function of 
the rods, f{9), or S = 2t[ j P2{cos{6))f{6)d6, where P2 is the second Legendre polynomial. 
For a completely aligned nematic S* = 1, whereas for an isotropic phase S* = 0. For rigid 



rods the predicted nematic order parameter at coexistence is S* = 0.79 



. The predictions 



from the Khokhlov-Semenov theory show quantitative agreement with the measured I-N 



transition for suspensions o: 



charged semiflexible virus fd, charged polymer xanthan, and 



neutral polymer PBLG J3|, |4 1 . 

Electrostatic interactions are incorporated into the Onsager model by rescaling the bare 
rod diameter D to a larger effective diameter Dgfr which depends on the ionic properties of 
the particle and the solution [l|, |9| . D^q is calculated from the second virial coefficient of 
Onsager 's free energy equation for charged rigid rods. In Fig. ^we plot Dcs as described by 
Stroobants et al. [9| as functions of ionic strength and rod surface charge. The non-linear 
Poisson-Boltzmann equation used in Stroobants description oi D^ was solved numerically 
using the approximations developed by Philip and Wooding [lO|. With increasing ionic 
strength D^s decreases approaching the bare rod diameter. Past experiments have shown 
that Deff accurately describes the ionic strength dependence of the I-N transition of fd virus 
suspensions J3[. For highly charged rods, the effect of surface charge on Z^efr is small as 
the non-linear nature of the Poisson-Boltzmann equation leads to counterion condensation 
near the colloid surface which renormalizes the bare surface charge to a lesser effective 
charge, which is nearly independent of the bare surface charge. In the nematic phase, the 
effective diameter increases due to an added effect called "twist" which is characterized by 
the parameter h = k,~^/Dcs, where k~^ is the Debye screening length. The effect of twist 
on Dcs, however, is predicted to be small for fd J3|, and we neglect it here. Studying the 
influence of ionic strength and surface charge on the I-N phase behavior tests if D^s can be 
accurately used to map charged rod phase behavior to hard-rod theories. 

Onsager's theory is based on an expansion of the free energy truncated at the second 
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FIG. 1: (Color Online) Effective diameter as a function of ionic strength and surface charge. With 
increasing ionic strength L>eff approaches the bare diameter of fd (M13) D = 6.6 nm. The effective 
diameter is plotted for surface charges of 10 e~/nm, 7 e~/nm and 4 e~/nm. These surface charge 
densities are the same as those measured for (a) fd at pH 8.2, (b) M13 at pH 8.2 or fd at pH 5.2, 
and (c) M13 at pH 5.2. At these surface charge densities DeflF is insensitive to variation in charge. 



virial level, so that only two-particle interactions are considered. This assumption has been 
shown to be accurate in the limit of very long rods, where L/D > 100 |ll|, or for very 
dilute suspensions. In our experimental system, however, decreasing the ionic strength 
rapidly decreases our effective aspect ratio to values far below the L/D = 100 limit. In 
order to accurately predict the phase behavior of rods with an effectively small aspect ratio, 
the theoretical free energy needs to incorporate third and higher virial coefficients. Scaled 
particle theory (SPT), which incorporates all higher virial coefficients in an approximate way 
. one theo^ wMch accomplish^ thifl. A scaled particle theo^ fo. ha.d „g.d .od. wa. 
originally developed by Cotte r | 12l| . More recently we have expanded this theory to include 
charge and semifiexibility J13l.ll4|. In conjunction with the Khokhlov-Semenov second virial 
theory we use this scaled particle theory to interpret our experimental results. 

In this paper we present experimental measurements of the isotropic-nematic phase tran- 
sition of semifiexible charged colloidal rods as a function of rod length, surface composition, 
solution pH, and solution ionic strength. We measure both the coexistence concentrations 
and the nematic order parameter and compare our results to both Onsager's theory, by 
way of Chen's numerical calculation!^], and scaled particle theory. For our model rods we 



use monodisperse suspensions of charged semiflexible rodlike fd virus, wild type M13 virus, 
and mutants of M13 virus which differ from the wild type only by their contour length. In 
solution, these particles exhibit isotropic, cholesteric (or chiral nematic), and smectic phases 
J3|, |l5|, |l6|, |l7|. Suspensions of fd have been previously shown to exhibit an I-N transition 
which agrees with theoretical predictions for semiflexible rods with an effective diameter Des 
[3|. M13 virus is structurally identical to fd virus, differing only in surface charge, making 
these two particles an ideal system for studying the influence of bare surface charge on the 
isotropic-nematic transition. Additionally, by comparing the I-N phase behavior of each 
of the M13 mutants, which except for length are structurally identical, and therefore by 
assumption have the same persistence length, we measure the influence of flexibility, defined 
as the ratio L/P on this transition. Though fd and M13 exhibit a cholesteric phase, the free 
energy difference between the cholesteric and the nematic phase is much smaller than the 
difference between the isotropic and nematic phases [l8|. This allows us to compare our 
results to theoretical predictions for the isotropic-nematic (I-N) transition. We refer to the 
cholesteric phase as the nematic phase henceforth. 

Motivation for these length and surface charge dependent measurements of the I-N tran- 
sition arose because new measurements of the nematic-smectic (N-S) transition in this same 

n 

system |13| exhibit measurable surface charge dependence and ionic strength dependence 
which can not be accounted for by treating the virus as a hard rod with a diameter D^s, 
in contrast to our previous measurements, which were limited in range of ionic strength 
[16[ . The new N-S measurements inspired a closer look at the ability of D^s to describe the 
effects of surface charge on the I-N transition. New measurements of the N-S transition as 
a function of length also indicate that semiflexibility has no measurable effect on the N-S 
transition for the limited range studied, which is as predicted, but which is in sharp contrast 
to the large predicted effect of flexibility on the I-N transition for the same range. The 
measurements presented here of the I-N transition as a function of charge and flexibility will 
contribute to the understanding of the relative importance of these variables in the evolution 
of the liquid crystalline ordering of charged semiflexible rodlike particles with concentration. 



II. MATERIALS AND METHODS 

Properties of fd and wild type M13 include length L =0.88 /im, diameter D = 6.6 nm, 
persistence length P =2.2 /im and molecular weight M = 1.64 x 10^ g/mol[20|. Each virus 
consists of approximately 2700 coat proteins helicoidally wrapped around single stranded 
DNA. The two viruses differ only by one amino acid per coat protein. In fd this amino 
acid is the negatively charged aspartate (aspi2), and in M13 it is the neutral asparagine 
(asnio)|2ll|. Thus at near neutral pH fd has one more negative charge per coat protein 
(3.4 ± 0.1 e~/protein) than M13 (2.3 ±0.1 e~/protein), which results in a net charge differ- 
ence of approximately 30% J23]. X-ray diffraction studies are unable to clearly discern any 
structural differences between M13 and fd J23|. The M13 length-mutants share the same 
properties as wild type M13, varying only in length and molecular weight, which scales 
linearly with length. The M13 mutant have lengths of 1.2/im, 0.64/im, and 0.39/im. Wild 
type Ml 3 J fd, and M13K07 (the 1.2/im mutant phage) were grown using standard tech- 



niques 



24| . The other two mutant phages were grown using the phagemid method, which 
produces bidisperse solutions of the phagemid and the M13K07 helper phage J2J]. We chose 
two plasmid DNA sequences, PGTN28 (4665bp) and LITMUS38 (2820bp) (New England 
Biolabs, Cambridge MA) to form our phagemids of length 0.64/im and 0.39/im, respectively. 
Sample polydispersity was checked using gel electrophoresis on the intact virus, and on the 
viral DNA. Except for the phagemid solutions, which contained approximately 20% by mass 
helper phage M13K07, the virus solutions were highly monodisperse as indicated by sharp 
electrophoresis bands. 

In a bidisperse system of long and short rods it is predicted that when isotropic and 
nemat ic p hases are in coexistence, the longer rods will strongly partition into the nematic 
phase J2^, |2a]. Using this fractionation effect we attempted to purify the bidisperse suspen- 
sions of the phagemid and M13K07 helper phage. We observed partitioning of the long rods 
into the nematic phase by DNA agarose gel electrophoresis (2-3 fold more long rods in the 
cholesteric phase than in the isotropic phase in qualitative agreement with Lekkerkerker et 
al. J2^), but were unable to successfully measure a difference in long rod concentrations 
in the isotropic phase after successive iterations of fractionation. The effect of fractiona- 
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TABLE I: Surface charge of fd and M13 at pH 8.2 and 5.2. (A) The charge of fd obtained by 
titration experiments [2j|. (B) M13 has one less negative amino acid per coat protein than fd, 
thus the surface charge of M13 can be approximated by subtracting one charge per protein subunit 
from the fd surface charge values. (C) Ratio of electrophoretic mobility (m), determined from Fig. 
121 of M13 to fd. (D) By multiplying the known fd charge by m, the linear surface charge density 
of M13 can be calculated. (E),(F) fd and M13 surface charge per unit length, respectively. 



tion on the coexistence concentrations was assayed by comparing the isotropic and nematic 
concentrations of coexisting samples (about 50% of each phase in one sample) with the 
highest concentrations for which the samples remained completely isotropic and the low- 
est concentrations for which the samples remained completely nematic, respectively. The 
only difference we observed was that the nematic concentration measured in coexistence 
with the isotropic phase was consistently about 5-10% lower than the nematic concentration 
measured when the sample was 100% nematic. The lower concentrations in the coexisting 
nematic phases are due to the partitioned long rods undergoing the I-N phase transition at 
lower mass concentrations. Because the effect of bidispersity is small, we report the phase 
behavior for the 0.39/im and 0.64 /xm rods at the limits of the coexistence region with the 
understanding that the samples contain about ~20% (by mass) 1.2 yum rods. 

All samples were dialyzed against a 20mM Tris-HCl buffer at pH 8.2 or 20mM Sodium 
Acetate buffer adjusted with Acetic Acid to pH 5.2. To vary ionic strength, NaCl was added 
to the buffering solution. The values for surface charge oi fd and M13 at pH 8.2 and pH 5.2 
are presented in Table HI The surface charge of fd was determined by titration experiments 



22J |. and the surface charge of M13 was calculated in two ways, both starting from the known 



fd surface charge. One way is to compare the molecular composition of fd and M13, and 
the second is to use the fact that because M13 and fd are identical except for their surface 
charge, their the electrophoretic mobilities are proportional to net surface charge ^2|. In 
Fig. we show using agarose gel electrophoresis of intact virus that fd migrates 200% faster 
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FIG. 2: Agarose gel electrophoresis of /d and M13 virus at pH 8.2 (gel 1) and pH 5.2 (gel 2). At pH 
5.2 the buffer was 40 mM Sodium Acetate, and at pH 8.2 the buffer was 40 niM Tris-Acetate-EDTA 
(TAE). Gels were run at ~1.0% agarose concentration and ~3.5 V/cm for 4 hours. Samples were 
placed in loading wells at a concentration of approximately 0.3 mg/ml. M13 and fd have the same 
length (L = 0.88/im) and diameter (D = 6.6nm), and differ only in surface charge. The ratio of 
electrophoretic migration distances between M13 and fd within each gel is therefore equal to the 
ratio of the surface charge. The electrophoresis bands for fd at pH 5.2 and M13 at pH 8.2 are not 
at the same migration distance, because the absolute migration distance is also a function of the 
buffer ions. 

than M13 at pH 5.2 and 150% faster at pH 8.2. Note in Table |l] we show that the surface 
charge of M13 at pH 8.2 is the same as the surface charge of fd at pH 5.2. 

All measurements were done at room temperature. The virus concentrations were mea- 
sured by absorption spectrophotometry with the optical density (A) of the virus being 
^269nm — 3-84 for a path length of 1 cm. The nematic order parameter was obtained by 
unwinding and aligning the cholesteric phase in a 2T permanent magnet (SAM-2 Humming- 
bird Instruments, Arlington, MA 02474) J2S| and measuring the sample birefringence. At 
2T, the magnetic field has a negligible effect on nematic ordering |2^|2I|. The nematic order 
parameters were calculated from the optical birefringence measurements obtained with a 
Berek compensator using the equation Arigat'S' = An, where Arigat is the saturation birefrin- 
gence. The value for Angat/p = 3.8 x 10^^[ml/mg], where p is the concentration of virus in 
[mg/ml] , as determined for fd via x-ray diffraction [l3[ . 
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FIG. 3: Isotropic-nematic coexistence concentrations as a function of M13 mutant contour length 
at 5 mM and 110 mM ionic strengths at pH 8.2. Open symbols represent the coexisting isotropic 
phase and solid symbols the nematic phase. Shaded areas are a guide to the eye indicating the 
coexistence regions. For rigid rods the coexistence concentrations pi oc 1/L at a constant ionic 
strength (constant I?cff)- Deviations from this relationship are most likely due to rod flexibility. 

III. RESULTS 

A. Effect of length and flexibility on the isotropic-nematic transition 



Figure El presents the length dependence of the I-N coexistence concentrations at high 
(110 niM) and low (5 mM) ionic strength. For rigid rods feeflfCj; the dimensionless concen- 
tration of the isotropic phase in coexistence with the nematic phase, is predicted to be a 
constant, feeflfQ = 3.29 J3l|, where 6eff = f-^^-Deg and Cj = piN^/M . In Cj, pi is the isotropic 
mass density, Na is Avogadro's number, and M is the molecular weight. Because the molec- 
ular weight is proportional to viral length, M = M^tL/ L^t, with M^t and L^^ equal to 
the molecular weight and length of wild type M13. Thus 6effQ = piLDcs{jLyjtN a/ M) = 
25pj[mg/ml]L[/xm]L'eff[/i m]. Therefore, for rigid rods, pi = const/L/Dcs, and at constant 
ionic strength (constant -Dcflf) Pi should be proportional to 1/L. However, we observe that at 
a given ionic strength, the slope of pj vs 1/L is not linear in Fig. IHl but instead increases with 
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FIG. 4: Dimensionless concentration of the isotropic phase in coexistence with the neniatic phase 
as a function of M13 mutant contour length for three ionic strengths at pH 8.2. The concentration 
is defined as ^effQ = j^Def[L?'Ni/V = 25/9i[mg/ml]L[/im]L'eflF[A'm]. Scale on the top of the graph 
identifies the flexibility in terms of L/P with P = 2.2//m. If the rods are rigid the phase behavior is 
predicted to be independent of length (Onsager) (dashed line). Semi- flexible rods show increasing 
6effCj with increasing flexibility as predicted by Khokhlov-Semenov theory calculated by Chen (solid 
line). Scaled particle theory at 100 mM ionic strength (SPTllO) and at 10 mM ionic strength 
(SPTIO) indicate that feeflfCj depends on L/Dcs- 



rod length, corresponding to an increase in feeffQ- This is shown more clearly in Fig. IH where 
feeflfCj is plotted as a function of length. The increase in feoffQ with length is in agreement with 
predictions for rods of increasing flexibility {L/P) , as shown by the theoretical curves from 
Khokhlov-Semenov theory and from SPT for semiflexible rods with a persistence length of 
P = 2.2/im. At high ionic strength ( I > 60 mM) we see good agreement with Khokhlov- 
Semenov theory calculated numerically by Chen (solid line) J6|. However with decreasing 
ionic strength, we measure an increase in the flexibility dependence of feeffQ- Subsequently, 
Khokhlov-Semenov theory only qualitatively describes the experimental results at low ionic 
strength. Agreement of the hard-rod Khokhlov-Semenov theory with our data is better at 
high ionic strength than at low ionic strength because the range of electrostatic interactions 
is weaker and L/Des is large, making the second virial approximation valid. 

To interpret the observed increase in flexibility dependence of the phase transition with 
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decreasing ionic strength, we turn to the scaled particle theory. The method for determin- 
ing the scaled particle theoretical coexistence concentrations and nematic order parameters 
is described elsewhere [lj| . In Fig. |3] we present the predicted SPT isotropic coexistence 
concentrations for rods with a diameter of 10.4 nm, (110 mM ionic strength), and 29.4 nm 
(10 mM ionic strength). At high ionic strength, SPT shows fair agreement with experi- 
mental results, and the theoretical curve for besCi is close to that predicted by Chen for 
the infinitely long rod limit. Additionally, we observe in Fig. |3]that SPT indeed predicts 
a small dependence of feesQ on L/Des, in contrast to the L/Des independent second virial 
theory. This suggests that effective aspect ratio of the rods, which decreases with ionic 
strength, has a small effect on the I-N transition concentration. However, the L/Deflf de- 
pendence predicted by SPT is opposite the trend experimentally observed; increasing Des, 
by lowering ionic strength, increases the measured besCi but lowers the scaled particle the- 
ory bcsCi. We argue that this discrepancy between scaled particle theory and experimental 
results at low ionic strength is due to the approximate treatment of electrostatics in Des, 
which is used not only as the theoretical hard rod diameter in SPT but also scales the ex- 
perimental coexistence concentrations from pi to bcsCi- -Dcflf is determined using from the 
second virial coefficient, and therefore is not necessarily accurate beyond that limit, ie. at 
low ionic strength. We note that the rescaled experimental coexistence concentrations, feesQ, 
are extremely sensitive to the value of Deff used to rescale the measured coexistence concen- 
trations. Pi. Differences in Des are translated linearly to changes in the experimental 6cffQ 
by bcsCi = 25pi[m.g/m\]L[fim\Des[p,'m]. However, the predicted effect of changing L/Dcs on 
feeflfCj, as shown by the SPT curves in Fig. ^ is much smaller than the measured change in 
besCi with ionic strength. Agreement between SPT and our experimental results improves 
if the effective diameter at low ionic strength is smaller than predicted at the second virial 
limit. 

The width of the coexistence region, (p„ — Pi)/pi, was also measured and is presented 
in Fig. El At low ionic strength, the coexistence width qualitatively follows the decrease 
expected for increasing flexibility shown by the solid line due to Chen J6|. For most rod 
lengths the value for the coexistence width is larger than predicted by both Khokhlov- 



11 



I 

cf 



0.6 



0.4 



0.2 



0.0 



l[mM] 



■ 


10 





60 


A 


110 




0.16 0.24 0.32 0.40 0.48 0.56 
UP 

FIG. 5: Width, (p„ — Pi)/ Pi, of the coexistence region as a function of rod flexibihty L/P. Results 
are plotted for three ionic strengths (10 mM, 60 mM, and 110 niM). Solid line is due to Chen for 
rods with P = 2.2/im 'G?|. Dotted and dashed lines are due to scaled particle theory for M13 rods 
(q=7e/nni) with a hard diameter Deff at 110 mM (SPTllO) and 10 mM (SPTIO) ionic strength, 
respectively. For rigid rods, the Onsager prediction for the I-N coexistence width is 0.29 |l|, la|- 
The width of the coexistence region should decrease with increasing flexibility. 



Semenov theory and by scaled particle theory. At short rod lengths this discrepancy is most 
likely due to the intrinsic bidispersity of the suspensions, which acts to widen the coexistence 
region [2^. A slow increase in the coexistence concentrations with time (possibly due to 
bacterial growth) pG] contributes to the large error bars, making comparison to predictions 
difficult. Above 10 mM ionic strength, where we see strong agreement between measurements 
of the coexistence concentrations and theoretical predictions, it is not apparent that there 
is any flexibility or ionic strength dependence in the width measurements. 

The nematic order parameter obtained from measurements of the birefringence of the 
magnetically unwound and aligned cholesteric phase in coexistence with the isotropic phase 
is presented in Fig. IHl We observe that at high ionic strengths, the nematic order param- 
eter decreases with increasing length (increasing flexibility) in qualitative agreement with 
Khokhlov-Semenov theory calculated by Chen p]. With decreasing ionic strength, however, 
the measured nematic order parameter increases, approaching Onsager's rigid-rod predic- 
tions, due to increasing the range of electrostatic interactions. This has also been observed 
for /o? virus suspensions |13|. Furthermore, at very low ionic strength (5 mM ionic strength) 
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the nematic order parameter becomes independent of rod length and equal to the predicted 
rigid rod value of 5 = 0.8. Scaled particle theory, as illustrated in Fig. IHl predicts that 
the nematic order parameter is largely independent of ionic strength. This suggests that 
the effective aspect ratio of the rods, which decreases with ionic strength, does not effect 
the nematic ordering. In addition, SPT agrees with the experimental measurements at high 
ionic strength better than Khokhlov-Semenov theory. 

Another possible explanation for an increase in nematic order parameter with decreas- 
ing ionic strength is electrostatic stiffening. If the interparticle interactions are dominated 
by electrostatics, the flexibility of the rods might be screened. This effective "electrostatic 
persistence length" P^i, which makes a charged polymer more rigid when in solution, is a 
dominant effect in determining the flexibility of charged flexible polymers with L/P ^ 1. 
However, for the semiflexible M13 and fd, Pei is predicted to be less than one percent larger 
than the bare persistence length ^3|. Additionally, the results for the coexistence con- 
centrations presented in Fig. |3] indicate that with decreasing ionic strength the measured 
coexistence concentrations deviate further from Onsager's rigid-rod predictions. Thus the 
measured coexistence concentrations and nematic order parameters exhibit contradictory 
trends, away from Onsager's rigid rod prediction versus towards Onsager's rigid rod predic- 
tion, respectively, with decreasing ionic strength. Therefore, electrostatic stiffening of the 
polymer cannot account for the observed high values of the order parameter at low ionic 
strength. Neither scaled particle theory, nor variation in the electrostatic persistence length 
satisfactorily explain the low ionic strength data. 

B. Effect of viral surface charge on the isotropic-nematic transition 

In this section we compare the phase behavior of M13 virus to that of fd virus as a 
function of surface charge and ionic strength. Recall that these particles have the same 
length L = 0.88/im and persistence length P = 2.2/im. In Fig. Owe present measurements of 
the isotropic coexistence concentrations as a function of viral surface charge at high and low 
ionic strength. The theoretical curve is from scaled particle theory for charged, semiflexible 
rods with L/P = 0.4. We only present the theoretical results from scaled particle theory 
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FIG. 6: (Color online) Nematic order parameter at coexistence as a function of rod length for four 
different ionic strengths. Solid black line represents the theoretical calculation by Chen ja] for the 
order parameter as a function of flexibility (L/P) indicated by the scale on the top of the graph. 
The dashed line is the theoretical nematic order parameter for rigid rods, S = 0.79[l|, |8!]. The 
scaled particle curves (dotted lines) are calculated as in J^] for virus rods at 110 mM (SPTllO) 
and 5 mM (SPT5) ionic strength. Theoretical curves were calculated for rods with a persistence 
length of 2.2/im. The measured order parameter decreases with increasing particle length at high 
ionic strength, but remains constant at low ionic strength. 
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FIG. 7: Coexisting isotropic phase concentrations pi as a function of particle surface charge for 
three ionic strengths, 10 mM, 60 mM, and 110 mM. Solid symbols are wild type M13 and open 
symbols are fd suspensions. Suspension pH is labeled above the graph for M13 and fd samples. Solid 
line is from scaled particle theory for semiflexible hard rods with a diameter Dog and L/P = 0.4. 
The charge dependence of the phase transition is well described by theory for ionic strengths of 60 
mM and 110 mM. 
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FIG. 8: (Color online) Width of the isotropic-neniatic coexistence region for wild type M13 and 
fd rods at three different surface charges as a function of ionic strength. Both M13 and fd have 
a length of L = 0.88/im. Solid symbols are at pH 5.2 and open symbols are at pH 8.2 for M13 
(circles) and fd (squares) suspensions. Solid line is from scaled particle theory for hard semiflexible 
rods with L/P = 0.4 and is independent of rod surface charge. The Onsager prediction for the 
I-N coexistence width in dimensionless units of be for hard rigid rods is (4.19 — 3.29)/3.29 = 0.29 
m, |g|. The coexistence width does not clearly show any charge dependence. 

in this section as this theory should more accurately describe the finite-length rod phase 
behavior than the second virial theory. In Fig. Owe confirm that the charge dependence of 
the I-N coexistence concentrations is accurately described by scaled particle theory at high 
ionic strengths. However, the efficacy of D^s as a means for incorporating all electrostatic 
interactions again diminishes at low ionic strength (I<60 mM), as seen previously in Fig. |3] 
and in Fig. El 

Fig. IHl presents the width of the coexistence region as a function of charge and ionic 
strength. The width of the coexistence region is independent of the surface charge of the 
rods and agrees (within large error bars) with scaled particle theory predictions. Both the 
measured coexistence concentrations and coexistence widths show that the effect of surface 
charge on the electrostatic interactions which drive the I-N phase transition are weak, which 
is consistent with the idea of charge renormalization incorporated into the calculations of 

Nematic order parameters obtained from measurements of the birefringence of magneti- 
cally unwound and aligned cholesteric samples of M13 at pH 8.2 are compared to previous 
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measurements of fd suspension nematic order parameters, measured via x-ray diffraction 
techniques, also at pH 8.2 |l3|, in Fig. El Recall that the nematic order parameter of 
fd is known to be proportional to the birefringence of the suspension by the relationship 
S = An/Ausat where Arigat = 3.8 x 10~^ml/mg |l3|. The order parameter of M13 was mea- 
sured at I-N coexistence as a function of ionic strength, and deep within the nematic phase 
for high (110 mM) and low (10 mM) ionic strength. Theoretical predictions from scaled 
particle theory for the nematic order parameter of hard semifiexible rods with L/P = 0.4 
are also shown in Fig. El The order parameters of M13 and fd were found to be equal as 
a function of ionic strength and concentration, indicating that the surface charge difference 
of 30% between the two particles does not effect nematic ordering. The insensitivity of the 
nematic order parameter to surface charge is consistent with the surface charge renormal- 
ization incorporated into D^s calculations (Fig. Q) [2|. The strong agreement of M13 and 
fd order parameters also indicates that these two different virus particles have the same 



3- 



Additionally, we again observe 



birefringence per particle, Arigat = 3.8 x 10~^ml/mg 

that the scaled particle theory fits the measured order parameter best for high ionic strength 

data. 

IV. CONCLUSION 

At high ionic strengths, where the range of electrostatic interactions are small and L/Dcs 
is large, the isotropic-nematic transition of the experimental system of charged semifiexible 
bacteriophages is well described by Khokhlov-Semenov theory for semifiexible charged rods. 
Increasing flexibility increases the coexistence concentrations besCi (Fig. H} and lowers the 
nematic order parameter (Fig. ^. In the region of high ionic strength, D^s accurately 
describes both the charge dependence and ionic strength dependence of the isotropic-nematic 
phase transition (Fig. |Z|). At low ionic strength, however, we find that the I-N coexistence 
concentrations and the nematic order parameter do not agree with theoretical predictions 
from either Onsager's second virial theory, or scaled particle theory. At low ionic strength, 
the flexibility dependence of the nematic order parameter is much weaker than expected 
(Fig. ini), but the flexibility dependence the coexistence concentrations is much stronger than 
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FIG. 9: Order parameter of the nematic phase (a) coexisting with the isotropic phase as a 
function of ionic strength (b)as a function of concentration at 110 mM ionic strength and (c) as 
a function of concentration at 10 mM ionic strength pH 8.2. Values for M13 were obtaine d by 
birefringence measurements and values for fd were obtained previously by x-ray diffraction |l3i |. 
Solid lines are scaled particle theory for semiflexible hard rods of diameter L^efj and L/P = 0.4. 
The order parameters for M13 agree with those measured for fd independent of concentration and 
ionic strength. 

expected (Fig. 0}. Because of these contradictory results we suggest that the disagreement 
between theoretical predictions and experimental data at low ionic strength is due to the 
approximate incorporation of the electrostatic interactions into the theoretical free energy 
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